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Spontaneous Charging and Crystallization of Water Droplets in Oil 
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We study the spontaneous charging and the crystallization of spherical micron-sized water-droplets 
dispersed in oil by numerically solving, within a Poisson-Boltzmann theory in the geometry of a 
spherical cell, for the density profiles of the cations and anions in the system. We take into account 
screening, ionic Born self-energy differences between oil and water, and partitioning of ions over the 
two media. We find that the surface charge density of the droplet as induced by the ion partitioning 
is significantly affected by the droplet curvature and by the finite density of the droplets. We 
also find that the salt concentration and the dielectric constant regime in which crystallization of 
the water droplets is predicted is enhanced substantially compared to results based on the planar 
oil-water interface, thereby improving quantitative agreement with recent experiments. 

PACS numbers: 68.05.-n, 64.70.D-, 82.70.Kj 
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I. INTRODUCTION 



It is "well-established that water and oil do not mix: 
droplets of water in oil (or droplets of oil in water) tend 
to coalesce such that the oil- water mixture coarsens until 
macroscopic phase separation of oil and water is achieved. 
It is also well-established that this coarsening process 
can be delayed or even prevented by additives such as 
surfactants or colloidal particles, which adsorb to the 
oil-water interface and thereby stabilize the droplets, ei- 
ther thermodynamically (such as in micro-emulsions) or 
kinetically (such as in Pickering emulsions) 1, 2]. Re- 
cently, however jexperimental observations by Leunissen 
and coworkers [3, |4j revealed stable micron-sized water 
droplets in somewhat polar oils without any additives. 
In fact, under appropriate conditions the oil-dispersed 
water droplets could even form a body-centered cubic 
(bec) crystalline phase with a lattice spacing of the or- 
der of 10 fim, and some of these crystals have been sta- 
ble for almost two years now without any observation of 
droplet coalescence [5j. The mechanism by which these 
water droplets are stabilized was argued to stem from 
an asymmetric partitioning of the (ever-present) mono- 
valent cations and anions over the oil and water phase; 
in the experiments of Refs. [3, Q with the oil cyclo- 
hexyl bromide the ions involved include H + and Br - , 
and the somewhat larger water affinity of the former com- 
pared to the latter should lead to positively charged wa- 
ter droplets. This suggested mechanism was confirmed 
in theoretical calculations [J, @ of monovalent ions in 
the vicinity of a planar oil-water interface, on the basis 
of Poisson-Boltzmann theory combined with ionic Born 
self energy in water and oil [7], [3 0- Although more 
advanced models could be invoked, e.g., involving more 
ionic correlations [1JJ, LLD, LL^i [l3| or a better account of 
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the ionic self-energy [Tdllfi Ha . [l7| , the relatively simple 
model used in Refs. [J, [|| showed surface charge densi- 
ties of the order of 10 — 100 elementary charges per /xm 2 
and hence 100 — 1000 charges for micron-sized water- 
droplets. With such a droplet charge the observed bec 
crystals of water-droplets in oil could be explained, at 
least qualitatively, and therefore we use this relatively 
simple Poisson-Boltzmann-Born model for further theo- 
retical explorations. 

The theory presented in Refs. 0, [fj considers a pla- 
nar oil-water interface separating two half spaces of oil 
and water. The advantage of this assumption lies in the 
fact that it allows for some analytic expressions for the 
surface charge density and the ionic contribution to the 
interfacial tension within nonlinear Poisson-Boltzmann 
theory [18J , which leads to an efficient scheme to analyze 
the parameter space. However, one could a priori expect 
quantitative shortcomings due to the assumed planar ge- 
ometry, e.g., because the typical experimental droplet ra- 
dius of about 1 /im is quite a bit smaller than the typ- 
ical screening length of about 10 /im in the oil phase, 
or because the typical lattice spacing in the crystal is of 
the same order as the screening length such that the net 
charge of the droplets could be affected by the nearby 
other droplets. In order to investigate these effects we 
extend in this paper the theory of Refs. 0, [fj] from the 
planar to the spherical geometry. This will be done in 
the context of a cell model [19l |. where a single spheri- 
cal droplet is considered in the center of a spherical cell 
with a finite volume representing the density of droplets. 
Whereas this geometry does no longer allow for analytic 
solutions of the Poisson-Boltzmann equation, the numer- 
ical solution is, however, fairly straightforward because of 
the radial symmetry. This geometry therefore enables us 
to study simultaneously the effects of droplet-curvature 
and droplet-density on the ionic double layers and on 
preferential adsorption and charging in the vicinity of the 
droplet surface. We will show that these effects, when 
compared to the planar limit results, give rise to a sig- 
nificantly larger crystallization regime for water-droplets 
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FIG. 1: A single Wigner-Seitz cell with radius _R centered on 
a droplet of radius a. To distinguish between a water-in-oil 
emulsion (WO,+) and an oil-in-water (OW,— ) emulsion, a 
sign convention has been introduced. 



in oil (due to a larger surface charge), and to a much 
smaller surface charge for oil-droplets in water. In fact, 
our numerical predictions for the crystallization regime 
are now quantitatively closer to the experimentally ob- 
served one, although there is still some deviation that we 
attribute to other shortcomings and oversimplifications 
of our microscopic model, e.g., the crude approximation 
of describing the ionic self-energies in oil and water by a 
simple Born energy. Our present results indicate, how- 
ever, that the essential physical mechanism of preferential 
ion partitioning can indeed explain the crystallization of 
water-droplets in oil as observed in Refs. [3, |j|. 



II. DENSITY FUNCTIONAL THEORY 
FOR SALINE EMULSIONS 

A. Wigner-Seitz cell approach 

We consider an emulsion of water-in-oil droplets (WO) 
of total volume V containing TV identical droplets with 
radius a. The volume of water in the system is defined 
as xV = 47riVa 3 /3, with x the volume fraction, hence 
(1 — x)V is the volume of oil. The theory describing 
emulsions of oil-droplets in water (OW) is analogous to 
that for WO systems outlined in this section. The emul- 
sion contains monovalent ions with ionic radii a±, which 
are typically 2 — 4 A. Using a cell model [19J , we reduce 
the iV-droplet problem to that of a single droplet in a 
spherically symmetric Wigner-Seitz cell, see Fig. [I] with 
radius R = ax^ 1 / 3 , such that V = AttNR 3 /3. 

In this single cell we consider a spherical oil-water in- 
terface located at r = a. The interface separates two 
bulk phases consisting of water (0 < r < a) and oil 
(a < r < R). Both, oil and water, are considered to be 
incompressible linear dielectrics, which means that the 
solvent background is characterized by the relative di- 
electric constant e w and e D , respectively, see Fig. [TJ The 



dielectric profile (relative to the dielectric constant of vac- 
uum e v ) is a step- function e(r) = e w if < r < a and 
e(r) = e if a < r < R. The ions are described by spheri- 
cally symmetric ionic density profiles p± (r) . The ions ex- 
perience Coulombic ion- ion interactions, which we treat 
in a mean-field fashion, and ion-medium interactions. 

The ion-medium interaction is taken into account via 
an external potential acting on the ions. Due to the di- 
electric properties of oil and water, the ions have different 
electrostatic self-energies in the two solvents. Using the 
Born approximation [9(, this self-energy of a cation (+) 
and an anion (— ) is given by E±(ei) = e 2 / '(2e v €ia±) , with 
e the elementary charge and i = w,o. This self-energy 
and the above dielectric profile e(r), allows us to rewrite 
the external potential acting on a cation and an anion 
as V±(r) — E±(e(r)) — E±(e w ), where we note that this 
potential is constructed to be zero in water. For realistic 
e « 4 — 20 the potential is of order (1 — 20)fcsT in oil, 
i.e., the ions prefer to be in the water. We also use the 
notation V±(r) = if < r < a and V±(r) = k B Tf± if 
a < r < R, with f± dimensionless and implicitly depen- 
dent on e w and e . Here fcgT is the thermal energy, ks 
Boltzmann's constant, and T the temperature. 



B. Poisson-Boltzmann equation 

Using the above external potential we employ the 
framework of density functional theory [20, [21j to cal- 
culate the equilibrium density profiles p±(r). The grand- 
potential functional fi[p±] for a single Wigner-Seitz cell 
can be written as 



r R 



log 



,fPi(r) 



+^QiHr,[p±}) + PVi(r) 



p, 
dr, 



- 1 



(1) 



with j3 = l/(fcflT) and the ionic valencies q± = ±1. The 
first line is the ideal-gas grand-potential functional. The 
second line describes the ion-ion Coulomb interaction 
in mean-field approximation and the ion-solvent interac- 
tions characterized by the external fields. The chemical 
potentials are represented in the form of an ion concen- 
tration p s , which is actually the ion concentration in a 
water reservoir in equilibrium with the emulsion. The 
electrostatic interactions between the ions in Eq. ([1]) are 
given in terms of the electrostatic potential functional 
ksT(f>(r, [p±])/e, which satisfies the Poisson equation 



e„e(r)V 2 <Kr, \p±\) 



with boundary conditions 



-4irf3e 2 ^2q l p l (r), 



lime w c/)'(r,[p±}) = lime o 0'(r, [p±]); 

rja rla 

lW(r,[p±]) = lW(r,[p±]) = 0, 

rj.0 rJH 



(2) 

(3) 

(4) 



where the prime denotes a derivative w.r.t. r. 

Minimizing the grand-potential functional leads to the 
Euler-Lagrange equation 5il/Sp±(r) = 0, which can be 
rewritten as Boltzmann distributions 



P±(r) = p s exp(-(3V ± (r) =p <l>(r, [p±]))- 



(5) 



In practice we implement the condition that p±(0) = p s 
for WO emulsions and p±(—R) = p s for OW emulsions, 
where we adhere to the sign convention explained in 
Fig. [TJ In the systems studied n w a ^> 1, with k -1 the 
Debye length in water, so that the water phase can in- 
deed be considered a bulk phase and hence acts as a salt 
reservoir with total ion concentration 2p s . Using Eq. ([5|) 
the Poisson equation reduces to 



(r) = n(r) 2 sinh(0(r) - 4> c {r)) (r ^ a), 



(0) 



where we have introduced n(r) = k w if < r < a and 
n{r) — k if a < r < R, with k 2 = 8irf3e 2 pi/(e v e i ) in 
medium i — o,w with p w — p s and p — p s exp(— [/ + + 
/_]/2) the bulk ion concentrations. The Donnan poten- 
tial <f> c (r) = (3[V-(r) - V+{r)]/2 follows from the local 
charge neutrality in the bulk liquids. Note that we have 
dropped the explicit [p±] dependence from 4>(r,[p±]) in 
Eq. ([6]) , since the above differential equation is not explic- 
itly P± dependent. Equation ^ together with its bound- 
ary conditions (see Eqs. $5§ and ([!])) has a unique solution 
and can be solved numerically on an r-grid by employ- 
ing standard numerical algorithms. Typically we require 
several thousand non-equidistant grid points, with a rel- 
atively small grid spacing close to r = a, tailored to the 
screening length in the oil and water phase, respectively. 



C. Ion-induced physical quantities 

Using the numerical solution for <p{r) (and hence for 
p±{r)) we can determine the charge and the excess inter- 
facial tension of the droplet, and the inter-droplet cou- 
pling parameter. The number of net unit charges induced 
by ion partitioning on the droplet is given by 



r a 



r 2 qiPi(r)dr. 



(J) 



Two charged water droplets, at separation r > 2a, are 
assumed to interact with each other through a screened 
Coulomb droplet-droplet interaction potential 



U(r) 



Z 2 e 2 exp(K (2a - r)) 



e v e (1 + K a) 2 r 
We introduce the coupling parameter 

T = (3U(p- 1/3 )(l + k + k 2 /2), 



(8) 



(9) 



with k = Kop- 1 / 3 and p -1 / 3 = (4 7 r/3) 1 / 3 i?, i.e., p = N/V 
the droplet density. It is empirically known from point- 
Yukawa {n a = 0) simulations 22] that crystallization 



occurs when T > 106. Even though the systems of 
present interest do have a finite 'hard' core, we can still 
apply this freezing criterion, because [/(2a) ^S> ksT and 
K a < 1 for most of our parameters, i.e., the repulsions 
are dominated by the screened-Coulomb part rather than 
the hard core. 

The excess interfacial tension 7 induced by ionic par- 
titioning is defined as 7 = (£![p±] — Vt[ph\) / {Ana 2 ) , where 
Ph{r) — p w if < r < a and Ph(r) — p if a < r < R, 
i.e., the difference between the grand-potential of the sys- 
tem and the grand-potential of a homogeneous reference 
system (per area). Using the above definition, 7 can be 
rewritten as 



1 






a ti 



Pi{r) - Ph(r) + T, c liPi{ r ) ( t>( r ) dr - 

(10) 

The term 'excess' refers to the fact that the total inter- 
facial tension reads 7tot = 7barc + 7, where 7baro is the 
bare oil- water interfacial tension, i.e., for an oil- water 
interface without the presence of ions. This bare interfa- 
cial tension is positive and typically of the order 1 — 10 
mN/m, whereas 7 turns out to be negative and of the 
order 10 - 100 nN/m. 



III. RESULTS AND DISCUSSION 
A. Choice of parameters 

From the experiments of Refs. [3, |j| we know that 
crystallization of water-in-oil droplets has been observed 
for the following parameters. Water has e w = 80.0 at 
T = 293 K. For water-droplets in CHB (cyclohexyl bro- 
mide) we find from an analysis of the published snap- 
shots that a = 1.1 ± 0.2 /zm and I = 9.6 ± 0.8 pm, with 
I the nearest neighbor distance in a bcc crystal Q. El- 
ementary geometry shows that the volume fraction of 
water is then given by x — 7rv3(a/Z) 3 , which yields 
x = (8.1 ± 1.6) • 10" 3 . For CHB e Q = 7.9 and electro- 
conductivity measurements give an indication of the bulk 
salt concentrations p w and p Q. According to Ref. [J] 
the major constituent ions are H + , OH - and Br - , with 
a H + = 2.8 A and a Br - = 3.3 A 3|. Likewise, for water- 
droplets in a CHB-decalin mixture (see Ref. [3|] for de- 
tails) we find a = 1.3 ± 0.2 pm, I = 16 ± 2 pm. and hence 
x = (2.9 ± 0.6) • 10~ 3 . This CHB-decalin mixture has 
e = 5.6 and k~ 1 > 3.6 pm [4]. Again the contributing 
ions are H + , OH - and Br - , but their respective concen- 
trations in CHB-decalin or water have a high degree of 
uncertainty. 

In our theoretical investigation we have chosen system 
parameters in the range indicated by the experiments of 
Ref. |3| . However, to capture the physics of the curvature 
effects present at spherical interfaces we do not fully take 
into account the complex chemistry described above. Our 
basis parameter set is R = 10 pm, a = 1 /im, p w = 10~ 3 
M, e w = 80, e = 5, a + = 3.6 A and a_ = 3.0 A. We vary 
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TABLE I: Physical quantities corresponding to Fig. [2] com- 
pared to those calculated analytically for a planar system. 
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FIG. 2: The double-layer near the interface of a spherical 
water- in-oil droplet, showing ion partitioning of the ions. Here 
a = 1 /im, p w — 10~ 3 M, e w = 80, e = 5, a+ — 3.6 A and 
a_ = 3.0 A, which gives re" 1 = 9.63 nm, re" 1 = 8.44 /an and 
p — 8.15 • 10" 11 M. The deviation from the homogeneous 
density profile is given as a function of the distance from the 
interface measured in screening lengths, for several values of 
the Wigner-Seitz cell radius R. The upper-left and lower-right 
quadrants correspond to p_ (r) profiles, whereas the lower-left 
and upper-right correspond to the p+(r) profiles. Note that 
the lines on the oil side of the interface terminate at (R—o)k , 
i.e., r = R. 



one or more of these parameters at a time and examine 
the effect on the physical quantities Z, 7, and I\ For this 
basis parameter set we have /+ = 14.8, /_ = 17.8 (using 
the Born approximation), x = 10~ 3 , p = 8.15 • 10 



11 



M, 



9.63 nm and n n 



.44 pm. We take two 



ionic species for simplicity and numerical convenience. 
The positive ionic radius corresponds to that of Na + and 
the negative to that of Cl~ [H, [H, HI ■ At this point 
the choice for the ion concentration in water seems a bit 
arbitrary, but we will show that it is in fact reasonable. 
The results for spherical interfaces have been calcu- 
lated using numerical techniques, whereas in the planar 
limit results can be determined analytically (la] . Note 
that the planar limit describes flat oil- water interfaces, 
therefore the planar system should correspond to a spher- 
ical Wigner-Seitz cell with a,R — a 3> k" 1 , i.e., the in- 
terface is locally flat on the scale of the Debye length of 
oil; we will denote this limit as a, R — ► 00. 



B. Preliminary analysis 

Figure. [2] shows the ionic density profiles for an a = 1 
/im water-in-oil droplet near the interface for several 
Wigner-Seitz cell radii R. Note the different scaling of 
all four axes, and that x increases from (R — ► 00), 
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to 0.001 (R = 10 A*m), to 0.008 (R = 5 /im), to 0.064 
(R = 2.5 lira). Ion partitioning causes the water phase 
to become negatively charged for a + > a-, whereas the 
oil phase picks up an equal but opposite charge. Note 
that local charge neutrality at r = R is violated for the 
finite Wigner-Seitz cells. This finite cell 'compresses' the 
double-layer inside the droplet w.r.t. that of the planar 
limit system, resulting in a droplet charge reduction, as 
can be seen from the shrinkage of the area enclosed by 
p+(r) and p-(r) at the water side in Fig. [U Of course, 
the same area-shrinkage occurs in the oil phase by the 
condition of global charge neutrality, however, this effect 
is not clearly visible. The value of the physical quanti- 
ties for the various systems in Fig. [2] is shown in Table U 
together with the values for the corresponding planar sys- 
tem. Note that we have introduced a = Z/(4ira 2 ), which 
is henceforth referred to as the 'surface charge density'' 
of the droplet. The total charge of a planar 'droplet' is 
infinite, because of its infinite surface area. The results 
in Table |T] show that the physical quantities Z, 7, and T 
are highly sensitive to the size R of the Wigner-Seitz cell 
in the experimentally relevant regime R *w 5 — 10 pm, 
e ~ 5, emphasizing the importance of curvature. We 
will examine this dependence more closely for extremely 
dilute emulsion, i.e., those emulsions which can be mod- 
elled by a droplet in an infinitely large Wigner-Seitz cell. 



IV. EXTREMELY DILUTE EMULSIONS 

A. Droplet size and curvature effects 

We now disentangle the effects of a finite droplet den- 
sity (finite R) and that of droplet curvature (finite) a by 
studying the curvature dependence in the extremely di- 
lute limit (R — > 00). In Fig.[3]the droplet's surface charge 
density a and the excess interfacial tension 7 as a func- 
tion of 1/a are shown in the extremely dilute limit. For 
the systems considered in Fig. [31 it follows that there is 
correspondence between the analytic planar results and 
the limit a — > 00, extrapolated from the spherical results, 
within a fractional uncertainty of w 10 -4 . Data was ob- 
tained for sufficiently large a to safely extend the lines 



through the point 1/a — 0. The surface charge density 
of the systems considered here (a > 100 nm) lies in the 
1—100 e/ jiva 2 range. We thus find that a = 1 /im droplets 
have charges between 10 — 1000 e. The excess interfacial 
tension ranges between 10 — 100 nN/m and is negative. 
The sign for this interfacial tension is a consequence of 
the external potential model we use here. We refer to 
Rcf. [26] for a more in-depth discussion of the excess in- 
terfacial tension of an interface separating two electrolyte 
solutions and the way in which it can be modelled. 

The curves in Fig. [3] show that there are two asymme- 
tries between OW and WO emulsions. Firstly, the devia- 
tion from the planar limit value is linear in 1/a in the case 
of WO systems (a > 0) and non-linear for OW systems 
(a < 0). This asymmetry can be explained entirely by 
the fact that K w a 3> 1 for the WO emulsions considered, 
whereas n a <C 1 for the OW emulsions. The devia- 
tion is linear for OW systems only in the small regime 
near 1/a = where K a > 1. For n w a « 1 (1/a ~ 100 
1/fim) the deviation becomes non- linear in WO emul- 
sions as well, however, this is well beyond the scale used 
in Fig. [3J Secondly, with decreasing droplet radius, \<j\ 
and I7I decrease for OW systems, but these quantities 
increase for WO emulsions. This effect can be attributed 
to the double-layer modification which occurs in spher- 
ical systems. It can be shown that the oil phase im- 
poses the total structure of the double-layer and hence 
determines the surface charge density and excess inter- 
facial tension. Curvature compresses the oil part of the 
double- layer in OW emulsions, resulting in a reduction of 
\a\ and (7], whereas the oil part of the WO double-layer 
gets stretched resulting in a corresponding increase of 
these quantities. Note that the regime in which the pla- 
nar limit approximation gives accurate results for a and 
7, i.e., a deviation of less than say 20%, is quite small. 
For an OW emulsion with e < 5 the droplet 's surface 
charge densityis negligible for experimentally reasonable 
droplet radii [3| . The planar value is an upper bound for 
the quantities a and 7 in extremely dilute OW emulsion. 
In WO emulsions, there may be an increase by 50% (or 
sometimes much more, see Table |TJ) for a = 1 fim. due to 
curvature effects. 



B. Curvature expansions 

It is known for fluids with all intrinsic (correlation) 
length scales smaller than all geometrical length scales 
that the deviation of any intensive quantity from the pla- 
nar value is a linear combination of mean and Gaussian 
curvature only 27]. In the spherical geometry consid- 
ered here this statement translates into a surface charge 
density as a function of the droplet radius a of the form 



a {a) — a p sign(a) + 



Cl 



C-2 



\n a\ \n a\< 



(11) 



with ci and c 2 coefficients and <r p the analytically known 
planar value for WO emulsions, provided k^ 1 ,^ 1 <C 
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FIG. 3: The surface charge density a (a) and the excess inter- 
facial tension 7 (b) for an extremely dilute system (R — > 00) 
with salt concentration in water p w = 10 -3 M, as a function 
of the water-droplet radius a at several relative dielectric con- 
stants e of oil. As specified in Fig. [l] positive values of a cor- 
respond to WO systems and negative values to OW systems. 
The inset shows a enlargement of the area close to the origin. 
The analytic planar limit values are indicated with dots. Note 
that the water-in-oil droplets have negative charge, whereas 
the oil-in- water droplets have positive charge. 



a, R. We use the sign convention introduced in Fig. [TJ 
which ensures that Eq. (fTTj) is valid for both OW and 
WO systems. A similar expression can be found for 7, 
by replacing a p with sign(a)7 p . Moreover, ci and C2 are 
positive and for typical system parameters of order unity. 
Equation (fTTj) proves useful as it allows us to describe the 
behavior of emulsion droplets for a range of droplet radii 
by determining the charge and excess interfacial tension 
for only two values of the droplet radius, fixing ci and C2 
in combination with the planar value. If the intrinsic and 
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FIG. 4: The surface charge density a (a) and the excess in- 
terfacial tension 7 (b) as a function of the relative dielec- 
tric constant of oil e for an extremely dilute emulsion, with 
p w - 10~ 3 M. Several droplet radii, a = 0.25, 0.5, and 1.0 
/jm, are compared with the analytic planar values (a — » 00). 
The charge of the droplets in OW emulsions is positive and 
in WO emulsions negative. The insets show the behavior of 
a and 7 for the larger domain 3 < e < 80, with both the 
planar and the a — 0.25 pm result (OW and WO) indicated. 
The OW line terminates due to numerical instabilities. 



the geometric length scales are not well separated, higher 
order terms in mean and Gaussian curvature can appear 
in Eq. ([TT]) [271 ] . However, these higher order terms turn 
out to be small in this study, since Eq. (fTTj) quantita- 
tively accounts for our numerical data in the droplet ra- 
dius regime |a| > 0.1 /zm even though k~ 1 > a in part of 
this regime. 



C. Varying e , p w and the droplet radius 

The effect of double-layer modification on a droplet in 
an extremely dilute system can also be evaluated when 
we vary e and p w in an experimentally reasonable range. 
In Fig. 01 a and 7 are given as a function of e , with 
R — » 00 and p w = 10~ 3 M, for several droplet radii a. 
From the insets, which show the full e -regime, we can 
see that the planar limit approximation is very accurate 
in a large e -range. However, it becomes apparent that in 
the experimentally relevant e -range there is a significant 
deviation from the planar value. In fact, for e < 7 we 
see that OW systems hardly experience any electrostatic 
effects (er « 0,7 w 0), whereas for WO emulsions such 
effects are much stronger than planar theory predicts. 
We find that a is of the order 1 — 100 e/pm 2 and 7 is of 
order 10 — 100 nN/m and negative. 

Again there is an asymmetry between OW and WO 
systems, which can be explained by the difference in De- 
bye length w.r.t. the droplet size. Note that the asym- 
metry between OW and WO becomes smaller when the 
dielectric constant of the oil increases and hence the De- 
bye length in oil decreases. Between the two limiting 
cases e — 1 and e = e w , for which both a and 7 are 
negligible in the former case and must vanish in the lat- 
ter case (in our model), we find an extremal value for a. 
This extremum is explained by an increase in p 0l but a 



decrease in f± for e D 



and vice versa for e n 



1. 



Note, however, that 7 seems to diverge for the smallest 
water-in-oil droplet that we consider here (a = 0.25 pm) 
in the limit of small e shown in the inset of Fig. [4Jd. 
A similar, non- vanishing behavior can be observed for a 
in the inset of Fig. 0^, however, we do not find an ap- 
parent divergence in a for the e -values considered here. 
This is in contradiction with the intuitive idea that our 
model should have negligible surface charge and excess 
surface tension in the e — > 1 limit. Our theoretical in- 
vestigation cannot exclude instabilities and uncertainties 
in the numerical algorithm used to solve for <f>{r) in the 
extreme e -regime. A more detailed evaluation of the 
limiting behavior for WO and OW emulsions with eo~l 
(e.g., water-droplets in air and air-bubbles in water re- 
spectively) with a spherical interface is beyond the scope 
of this work. 

In Fig. [S] we show a and 7 as a function of p w in an 
extremely dilute system with e = 7.5 for several droplet 
radii a. One can see that the behavior of the physi- 
cal quantities for spherical interfaces w.r.t. their planar 
counterparts is analogous to that found in Fig. 0] This 
analogy can be easily explained by the way in which the 
Debye lengths are modified when changing either e or 
p w . Note the values of a and 7 are in the range 1 — 100 
e/pm 2 and 10 — 100 nN/m, respectively. Our results for 
extremely dilute emulsions, Figs. [21 LH and[5l thus show 
that for reasonable choices of the system's parameters 
I7I is of the order 10 — 100 nN/m. This excess interfacial 
tension therefore does not contribute significantly to the 
bare interfacial tension of an oil-water interface, which is 
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FIG. 5: The surface charge density a (a) and the excess inter- 
facial tension 7 (b) as a function of the bulk ion concentration 
in water p w for an extremely dilute emulsion with e — 7.5. 
Several droplet radii, a = 0.25, 0.5, and 1.0 pm, are consid- 
ered and can be compared to the planar value. The charge of 
the droplets in OW emulsions is positive and in WO emulsions 
negative. Some lines terminate due to numerical instabilities. 



of the order 1 — 10 mN/m. 



V. CRYSTALLIZATION AT FINITE DROPLET 
VOLUME FRACTION 

In this section we consider a Wigner-Seitz cell with 
R = 10 /im, yielding a bcc nearest neighbor distance 
I = (\/3tt) 1 / 3 R = 17.6 pm, which for e„« 5 is within 
the regime for the experiments of Refs. [3, |4(. Only WO 
emulsions are examined, as we found no crystallization 
for OW emulsions, i.e., T <C 106 (see Eq. ©) for any 



reasonable choices of the OW system parameters. Our 
model thus predicts that crystallization of oil-in-water 
droplets does not occur or is extremely unlikely. Figure [6] 
indicates for WO systems the T — 106 isoline, which is 
in fact the WO droplet freezing line, in the (e ,a;)-plane 
(a) and the (yOu,,a;)-plane (b) for our basis parameter set. 
The choice for the e -range in Fig. [5^ is inspired by the 
range of dielectric constants for which droplet crystal for- 
mation has been observed [3| . The p^-range in Fig. [6)3 
is physically reasonable, inspired by the isoline minimum 
found in planar analysis, and limited by the stability of 
our numerical algorithm to solve for <f>(r) in the spherical 
geometry. 

We see that there is an isoline minimum at eo ~ 5.5 
(p w = 10" 3 M) and at p w « 3.0 • 10~ 3 M (e = 5) 
with x ~ 0.025 for the spherical results. This value is 
substantially smaller than that of the planar minimum 
(x ~ 0.065), however it is still significantly larger than 
the experimentally found water content of emulsions in 
which water-in-oil crystals were observed. This was to 
be expected if one considers the uncertainty of some of 
the parameters used, particularly the ionic contents of 
the emulsions and the corresponding ionic self-energies. 
The location of the minimum also gives an a posteriori 
justification of our choice to use p w = 10~ 3 M for our 
basis parameter set. Note the regime in which crystal- 
lization can occur according to spherical theory is greatly 
extended with respect to that found using planar theory. 
Yet there are parameter choices for which this regime 
is reduced, also see Fig. [6] Therefore, we must con- 
clude that the effects of curvature on the crystallization 
of spherical water-in-oil droplets are non-trivial and in- 
volve competing processes: at extreme dilution the sur- 
face charge density is higher than the planar value, but at 
finite concentration it becomes smaller, also see Table fl] 
In accordance with the rough location of the minimum 
found in Fig. [SJ we will use e = 5 and p w = 10~ 3 M to 
examine the effects of the self-energy difference between 
water and oil of the respective ions, keeping R = 10 pxsx 
(I = 17.6 pm) fixed. 

In Fig. [7£i the T = 106 isolines are indicated as a func- 
tion of f± for several p w and x = 10~ 3 , together with the 
convex envelope of the isolines within the p w = 10~ ' — 10 
M range for both spherical and planar interfaces. The 
choice of this p w envelope range is inspired by well-known 
numbers for the ion concentration in water, which is 
bounded from below by that of pure water with a pH of 
7 caused by self-dissociation of water molecules and from 
above by that of a saturated solution of order 10 M. Note 
that crystallization of water-droplets in oil can occur in 
the region enclosed by the isolines or envelopes, respec- 
tively, and the /+-axis, i.e., in this region T > 106. The 
spherical interface envelope in Fig. [7^, was determined by 
calculating the T = 106 isolines for p w = 10~ 6 - 10~ 2 M. 
All of these isolines turn out to be shifted in the same 
way with respect to their planar counterparts as is shown 
explicitly for p w — 10 -3 M. Hence, we have assumed that 
this behavior can be extrapolated to p w = 10 M. Note 
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FIG. 6: The freezing line, F = 106, as a function of the relative 
dielectric constant of oil e and the water volume fraction x 
for R = 10 pra and p w — 10 -3 M (a) and as a function of p w 
and a; for R = 10 /im and e = 5 (b). The planar T = 106 
isoline is also indicated and several data points which lie on 
the spherical interface freezing line are included. Above the 
isolines F > 106 and below T < 106. Note the shift in the 
freezing line minima and the increase in the crystallization 
zones. 



again that the effect of curvature is to increase the range 
over which crystallization can occur. 

Fig. [TJd shows these convex envelopes, in the planar 
limit approximation, as a function of f± for several x and 
e . In agreement with our findings in Fig. [SI we recover 
that crystallization occurs more easily for larger x, as one 
would expect, since larger droplets have a higher charge 
and are closer together when R is fixed. Note that in the 
2.5 - 10 e -range there is no significant shift in the con- 
vex envelopes and their corresponding r = 106 isolines. 
The spherical interface results have not been considered 



FIG. 7: (Color online) Coupling parameter isolines V — 106 
for water-in-oil droplets in the (/+,/_) plane (see text), in 
(a) for several salt concentrations p w (dashed) and the p w - 
envelopes (full curves) based on both curved and planar ge- 
ometries at composition x = 10~ 3 and oil dielectric constant 
e — 5, in (b) only the p^-envelopes at several x and e based 
on the planar geometry. In all cases the droplet density is 
kept fixed such that R = 10 pm. The diagonal (/+ = /_) 
is a reflection symmetry axis for all curves. Crystallization is 
predicted between the /+ axis and the curves (and between 
their reflected images of course) . 



here, because of the time consuming character of these 
calculations. However, one can expect an increase in the 
crystallization zone for these envelopes similar to that of 
Fig. [Hi. 



VI. CONCLUSION AND OUTLOOK 

We have presented calculations for anions and cations 
near a spherical water-oil interface, taking into account 
ionic self-energies and screening, to describe the sponta- 
neous charging of water droplets in oil. This theory was 
applied to emulsions of oil and water which contain ions. 

In the extremely dilute droplet limit, the effects of cur- 
vature on the charge and excess interfacial tension in- 
duced by ion partitioning of the anions and cations were 
compared to results obtained for a planar interface. It 
turns out that the planar limit approximation used in 
Ref. p can be applied with a high degree of accuracy 
for many system parameters. However, in the range 
of the experiments of Refs. [3, [j], we have shown that 
spherical and planar results differ significantly. Water-in- 
oil droplets have a substantially higher and oil-in-water 
droplets a substantially lower charge/excess interfacial 
tension than one would expect on the basis of planar cal- 
culations. In accordance with Ref. [23] we found that the 
value of physical quantities in a spherical system can to 
an extent be approximated using a polynomial expansion 
in l/|«; a| around the planar value. 

For finite volume fractions of water in oil we have in- 
vestigated the crystallization of water-in-oil droplets us- 
ing the dimensionless coupling parameter T of a point- 
Yukawa system to predict crystallization [22| . The range 
in parameter space in which crystal formation can occur 
is greatly extended by using spherical values w.r.t. the 
planar result, mainly because of the larger surface charge 
densities in the spherical case. We expect that the the- 
ory we have presented captures the physics of the experi- 



ments performed by Leunissen et al. in Refs. 0,0]. How- 
ever, quantitative comparison between theory and exper- 
iment is not possible at this time. Not only additional 
theoretical effort is required, e.g., including more realistic 
self-energies [17J or a wider class thereof [lj, [l5|, [l6| , but 
the complex chemistry in these oil- water emulsions needs 
to be further scrutinized and elucidated experimentally 
to facilitate such a comparison. 

Another extension of the present theory for ions in the 
vicinity of a curved oil-water interface is the addition 
of charged colloidal particles. It was shown experimen- 
tally [3, 0] and theoretically in the planar geometry 0, [y] 
that the phenomenology due to the presence of charged 
colloids is extremely rich, e.g. involving strong colloidal 
adsorption at the droplet surface with adjacent huge 
colloid-free zones. One would expect curvature effects 
in these systems as well. Moreover, the much smaller oil- 
in-water droplets in the 10— 100 nm range as observed in 
the Pickering emulsions of Refs. [28|,[29| certainly warrant 
a theoretical treatment that takes the finite curvature of 
the droplets into account. Studies along these lines are 
in progress. 
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